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A systematic method is proposed for the determination of the statistical properties of a field 
consisting of a coherent structure interacting with turbulent linear waves. The explicit expression 
of the generating functional of the correlations is obtained, performing the functional integration on 
a neighbourhood in the function space around the soliton. The results show that the non-gaussian 
fluctuations observed in the plasma edge can be explained by the intermittent formation of nonlinear 
coherent structures. 



I. INTRODUCTION 

In a recent work it has been proposed a systematic analytical method for the investigation of the statistical 
properties of a coherent structure interacting with turbulent field. The method is here developed in detail and new 
possible applications or developments arise. 

The nonlinearity of the dynamical equations of fluids and plasma is the determining factor in the behaviour of these 
systems. The current manifestation is the generation, from almost all initial conditions, of turbulent states, with an 
irregular aspect of fluctuations implying a wide range of space and time scales. The fluctuations seem to be random 
and a statistical characterization of the fluctuating fields is appropriate. However it is known both from theory and 
experiment that the same fields can have, in particular situations, stable and regular forms which can be identified 
as coherent structures, for example solitons and vortices. For most general conditions one should expect that these 
aspects are both present and this requires to study the mixed state consisting of coherent structures and homogeneous 
turbulence. 

Numerical simulations of magnetohydrodynamics show that in general cases a coherent structure emerges in a 
turbulent plasma, it moves while deforming due to the interactions with the random fields around it and eventually 
is destroyed. In plasma turbulence a coherent structure is build up by the inverse spectral cascade or by merging 
and coalescence of small-scale structures ||], H), The nonlinearity of the equations for the drift waves in a 
non-uniform, magnetized plasma permits the formation of solitary waves in addition to the usual small-amplitude 
dispersive modes. The convective nonlinearity (of the Poisson bracket type), can lead to low-frequency convective 
structures in magnetized plasma |^ , ^ , |0] , g] , The structures are not solitons in the strict sense but are very 
robust. It is even possible that the state of plasma turbulence can be represented as a superposition of coherent vortex 
structures (generated by a self-organization process) and weakly correlated turbulent fluctuations. 

Naturally, the coherent structures influences the statistical properties of the fields (the correlations), in particular 
the spectrum. In this context it is usual to say that the deviation of the correlations of the fluctuating fields from the 
gaussian statistics is associated with the presence of the coherent structures and it is named intermittency. Numerical 
simulations ||l0|| of the 2-dimensional Navier-Stokes fluid turbulence have shown coherent structures evolving from 
random initial conditions and in general energy spectra steeper than have been atributed to intermittency (patchy, 
spatial intermittent paterns). These coherent structures are long-lived and disappear only by coalescence, the latter 
being manifested as spatial intermittency. In these studies it was underlined that the coherent structures have effects 
which cannot be predicted by closure methods applied to mode-coupling hyerarchies of equations. 

The difficulty of the analytical description consists in the absence from theory of well established technical methods 
to investigate the plasma turbulence in the presence of coherent structures. While for the instability-induced turbu- 
lence (via nonlinear mode-coupling) systematic renormalization procedures have been developed, the problem of the 
simultaneous presence of coherent structures and drift turbulence has not received a comparable detailed description. 

In the recently proposed method , the starting point is the observation that the coherent structure and the drift 
waves, although very different in form, are similar from a particular point of view: the first realizes the extremum 
and the later is very close to the extremum of the action functional that describes the evolution of the plasma. 
The analytical framework is developed such as to exploit this feature and is based on results from well established 
theories: the functional statistical study of the properties of the classical stochastic dynamical systems (in the Martin- 
Siggia-Rose approach); the perturbed Inverse Scattering Transform method, allowing to calculate the field of perturbed 
nonlinear coherent structures; the semi-classical approximation in the study of the quantum particle motion in multiple 
minima potentials. 
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The dilute gas of plasma solitons has been studied by Meiss and Horton |13|] who assumed a probability density 
function of the amplitudes characteristic of the Gibbs ensemble. We analyse the same nonlinear equation but take 
into account the drift wave turbulence. 

A brief discussion on the closure methods developed in the study of drift wave turbulence provides us the argumen- 
tation for the need of a different approach (Section 2). The Section 3 contains a description of the general lines of 
the method proposed. A more technical presentation of the calculation is given in the subsection 3.2. The particular 
case of the drift wave equation is developed in detail in the Section 4 and in the Section 5 the explicit expression of 
the generating functional is used for the calculation of the correlation functions. The results and the conclusions are 
presented in the last Section. Some details of calculations are given in the Appendix. 



II. THE NONLINEAR DYNAMICAL EQUATIONS 

We consider the plasma confined in a strong magnetic field and the drift wave electric potential in the transversal 
plane (x, y) where y corresponds to the poloidal direction and x to the radial one in a tokamak. We shall work with 
the radially symmetric Flierl-Petviashvili soliton equation |^ studied in Ref. : 

(l-p.Vj- + ..--«.^-=0 (1) 

where ps — Cg/i^i, Cg = {Te/rrii)^^^ and the potential is scaled a,s ip — -j^ |f . Here L„ and Lt are respectively the 
gradient lengths of the density and temperature. The velocity is the diamagnetic velocity Vd — ^jf^- The condition 
for the validity of this equation is: {k^Ps) {kps)^ <C Ve-ff-^ where rj^, — 4^- 



The exact solution of the equation is 



'Ps{y,t;yo,u) = -3 [ 1 sec/i^ 7^ {y - yo - ut) 



(2) 



where the velocity is restricted to the intervals u > Vd or u < 0. The function is represented in Fig.l. In the Ref. 

the radial extension of the solution is estimated as: (Ax)^ ^ PsLn- In our work we shall assume that u is very 
close to Vd , Vd (i.e. the solitons have small amplitudes). 

The nonlinear equations for the drift waves are known to generate as solutions irregular turbulent fields but also 
exact coherent structures of the type (||), depending on the initial conditions. Typical statistical quantities are the 
correlations, like: (ip {x, y, t) ip {x' , y' , t')) ~ |x — x'|^ \t — t'\^ , where for the homogeneous turbulence the exponents 
C and z are calculated by the theory of renormalization or by spectral balance equations, using closure methods ( 
pOf). Various closure methods have been developed as perturbations around gaussianity and they are valid for small 
deviation from the gaussian statistics . We see intuitively that this approach cannot be extended to the description 
of the coherent structures. This can also be seen in more analytical terms. A quantity which is unavoidable in the 
calculation of the correlations is the average of the exponential of a functional of the fluctuating field, consider simply 
(exp (p)) (for example in the inverse of the Vlasov operator, using the Fourier transformation, the potential appears 
in the formal expression of the trajectory, i.e. at the exponent). This quantity can be written schematically as p7|: 



(exp(v5)) = exp 



(3) 



where ((<^")) represents the cumulant of order n (i.e. the irreducible part of the correlation, after substracting the 
combinations of the lower order cumulants). For a gaussian statistics the first two cumulants are different of zero 
(n = 1 : average and n = 2 : dispersion), all others are zero. Non- vanishing of the higher order cumulants is the 
signature of non-gaussian statistics. In the perturbative renormalization we assume slight deviation from gaussianity, 
i.e. small absolute values of the next order cumulants (e.g. the kurtosis must be close to 3, the gaussian value) and 
vanishing of the higher order cumulants. This assumption is obviously invalid in the case of coherent structures. The 
field of a coherent structure has long range, persistent correlations imposed by its regular geometry, which naturally 
requires non- vanishing very large order cumulants (i.e. many terms in the sum at the exponent in Eq.(^)) and excludes 
any perturbative expansion. 

In particular, the closure of the nonlinear equation for the two-point correlation (based on the retaining the directly 
interacting triplet) can account for the small scale correlations related to the space-dependent relative diffusion, i.e. 
the clump effect ( ||2^), but the spectrum obtained in this framework cannot account for the possible 

existence of the coherent structures. This clearly suggests that we must find a different approach. 
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III. COHERENT STRUCTURES IN A TURBULENT BACKGROUND 



A. The outline of the method 

We present the basic lines of an approach which can provide a statistical description of the coherent structure in 
a turbulent background. The physical origin of this approach is the observation that the non-linear equation whose 
solution is the coherent structure (the vortex soliton) also has classical drift waves as solutions, in the case of very 
weak nonlinearity. In a certain sense (which will become more clear further on), the vortex soliton and the drift waves 
belong to the same family of dynamical configurations of the plasma. Our approach, which is designed to put in 
evidence and to exploit this property, consists of the following steps. 

We start by constructing the action functional of the system. The dynamical equation is the Euler-Lagrange 
equation derived from the condition of extremum of this functional and the exact solution is the vortex soliton (||) . 

By using the exponential of the action we construct the generating functional of the irreducible correlations of Lp. 
This functional contains all the information on the coherent structure and the drift turbulence. The correlations 
are obtained via functional differentiations. This requires the formal introduction of a perturbation of the system, 
through the interaction with an external current. Throughout the work, this perturbation will be considered a small 
quantity and finally it will be taken zero. 

The generating functional is by definition a functional integral over all possible configurations of the system and 
this integral must be calculated explicitely. The simplest thing to do is to determine the configuration of the system 
(with space and time dependence) which extremises the action, by equating the first functional variation of the action 
with zero and solving this equation: this will give the vortex soliton (modified due to the small interaction term). 
Then one should replace this solution in the expression of the action. This is the lowest approximation and it does 
not contain anything related to the drift wave trubulence. 

At this point we can benefit of the particular physics of the drift waves. The vortex soliton is the exact solution of 
the fully nonlinear equation and is a localized potential perturbation with regular, cylindrical symmetric form. The 
linear drift waves are harmonic potential perturbations propagating with constant velocity (the diamagnetic velocity 
in the case of the drift poloidal propagation in tokamak). Although the drift waves have very different geometry they 
are solutions of the same equation as the vortex, but for negligible magnitude of the nonlinear term. The drift waves 
do not exactly realize the extremum of the action functional, but obtain an action very close to this extremum. This 
means that the drift waves and the vortex soliton are close in the function space in the sense of the measure defined 
by the exponential of the action. In other terms the drift waves are in a functional neighbourhood of the vortex (for 
this measure). This suggests to perform the functional integral with better approximation, which means to perform 
the integration over a functional neighbourhood of the vortex solution. This will automatically include the drift waves 
in the generating functional of correlation which so will contain information on both the coherent structure and the 
drift waves. The function space neighbourhood over which the functional integration is extended is limited by the 
measure (exponential of the action) which severly penalizes all configurations of the system which are far from the 
solution realizing the extremum (i.e. the vortex soliton). As in any stationary phase method there are oscillations 
which strongly suppress the contribution of the configurations which are far (in the sense of the measure) from the 
soliton. In practice we shall expand the action in a functional Taylor series around the soliton solution and keep the 
term with the second functional derivative. 

In this perspective the drift waves appear as fluctuations around the soliton solution. This is compatible with the 
numerical simulations which show that the vortices are accompagned by a tail of drift waves. During the interaction 
of the vortices linear drift waves are "radiated" ||2^. On the other hand, the analytical treatment of the perturbed 
vortex solution by the perturbed Inverse Scattering Transform shows similar tail of perturbed field, following the 
soliton. This strengthens our argument that integrating close to the vortex means to include the drift waves in the 
generating functional. 

The functional integral can be performed exactly and we determine the generating functional of the potential 
correlations. We shall calculate the two-point correlation by performing double functional derivative at the external 
current. 



B. Expansion around a soHton 

1. The action and the generating functional of the correlations 

The analytical framework is similar to the model of quantum fluctuations around the instanton solution in the 
semi-classical calculation of the transition amplitude for the particle in a two- well potential (see reference Is^). Let 
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us write formally the equation for a nonlinear plasma waves as 



Oip = 



(4) 



where the field ip {x, y, t) represents the "field" (coherent structure + drift waves) and the operator is the nonlinear 
operator of the equation (|l|) . This equation should be derived from the condition of extremum of an action functional 
which must reflect the statistical nature of our problem. The field ip obeys a purely deterministic equation, but the 
randomness of the initial conditions generates a statistical ensemble of realizations of the system evolutions (space- 
time configurations). We shall follow the Martin-Siggia-Rose method of constructing the action functional but in 
the path-integral formalism, for which we give in the following a very short description ( [Q). First, we 

consider a formal extension from the statistical ensemble of realizations of the system's space-time configurations to 
a larger space of functions which may include even non-physical configurations. Every function is discretized in space 
and time, so it will be represented as a collection of varables (pi , each attached to the corresponding space-time point 
i. In this space of functions, the selection of the configurations which correspond to the physical ones (solutions of 
the equation of motion) is performed through the identification with Dirac delta-functions, in every space-time point 



(5) 



and integration over all possible functions ip, i.e. over the ensemble of independent variables ifi. Using the Fourier 
representation for every S function we get 



Y[dip, / Y[dx, exp \ixi0ip{xi,yi,ti) 



(6) 



Going to the continuum limit, a new function appears, x i^: Vi t) which is similar to the Fourier conjugate of ip. The 
generating functional of the correlation functions is 

D [ip (x,t)] D [x (x,i)] exp |z J d^'dt'x (x', t') 0^ (x', t')| (7) 

where the functional measures have been introduced and x = {x, y). 

The random initial conditions ipo (y) can be included by a Dirac 5 functional: 5 {ip (to, y) — ipo (y)). As explained in 
[01, instead of this exact treatment (accessible only numerically) we exploit the particularity of our approach, i.e. the 
connection between the functional integration and the delimitation of the statistical ensemble: the way we perform 
the functional integration is an implicit choice of the statistical ensemble. We choose to build implicitely the statistical 
ensemble, collecting all configurations which have the same type of deformations (given in our formulas by xj)- -'^H 
these configurations belong to the neighbourhood of the extremum in function space and we take them into account, 
by performing the integration over this space. In doing so we assume that the ensemble of perturbed configurations 
induced by an "external" excitation (J below) of the system is the same as the statistical ensemble of the system's 
configurations evolving from random initial conditions. 

We must add to the expression in the integrand at the exponential a linear combination related to the interaction 
of the fields tp and x with external currents J^p and J^: 



Sj = dx'dt 



D[ip{x,t)] D[xix,t)]e^p{iSj} 

X(x',t') d^ix',t') + J^ip + J^x 

It is now possible to obtain correlations by functional differentiation, for example 

1 Z 

{•^{x2,y2,t2)'~p{xi,yi,ti)) 



(8) 



(9) 



j=o 



Zj SJ^ ix2,y2,t2)SJ^ {xi,yi,ti) 
For the explicit calculation of the generating functional we need the functions ip and x which extremize the action 

6Sj 
Sip 
SSj 



(10) 
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2. Schema of calculation of the generating functional 



In the absence of the current J the equations (|T^) have as solutions for ip the nonhnear sohtons (vortices) |3^] , . 
More generally, the basic solution of the KdV equation (on which the Flierl-Petviashvili equation can be mapped) is 
the periodic cnoidal function which becomes, when the modulus of the elliptic function is close to 1, the soliton. When 
the distance between the centres of the solitons is much larger than their spatial extension (dilute gas) the general 
solution can be written as a superposition of individual solitons, with different velocities and different positions [p^ . 
For simplicity we shall consider in this work a single vortex soliton and in the last Section we shall comment on the 
extension of the method to many solitons. 

The position of the centre of the soliton rises the difficult problem of the zero modes [ ^2| . Except for a brief comment 
about the relation of the zero modes with the gaussian functional integration (see below) , we shall avoid this problem 
and postpone the discussion of this topic to a future work. 

In the presence of the external current J, the equations resulting from the extremization of the action Sj become 
inhomogeneous, and the solutions are perturbed solitons. This point is technically non-trivial and we shall use the 
results obtained by Karpman ||2^ who considered the Inverse Scattering Transform method applied to the perturbed 
soliton equation. We find the approximate solution ipjs and xjs of the inhomogeneous equations (i.e. including the 
external current J). The result depends on the currents J, and this will permit us to perform functional differentiations 
in order to calculate the correlation, as shown in Eq.(|^). As a first step in obtaining the explicit form of Zj, the 
perturbed soliton solutions depending on J must be introduced in the expression of the action Sj . After that we 
perform the expansion of the functions Lp and x around the coherent solution. 



V = Vjs + Sip 
X = XJs + 5x 



(11) 



This gives 



Z.j = e^^{iSjs) J D[Sp]D[Sx] 

X exp < / dx'dt' 6x(x',t') , , , 
1 J \ SpSx 



Sp (x',i') 



or 



Zj - exp {zSj,) ^ (27r)"/^ ( det 



dp6x 



-1/2 



(12) 



VJs.XJs 

since the integral is gaussian [p4[. The determinant is calculated using the eigenvalues 



SpSx 



and 



-0„ (X, t) = Xnlpn (X, t) 



(13) 



det 



SpSx 



(14) 



Since the action is invariant to the arbitrary position of the centre of the soliton there are directions in the function 
space where the fluctuations are not bounded and in particular are not Gaussian. This requires the introduction 
of a set of collective coordinates and after a change of variables the functional integrations along those particular 
directions are replaced by usual integrations over the colective variables, with inclusion of Jacobian factors. The 
zero eigenvalues of the determinant (corresponding to the zero modes) are excluded in this way. We shall avoid this 
complicated problem and assume a given position for the centre of the vortex. 
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IV. APPLICATION TO THE VORTEX SOLUTION OF THE NONLINEAR DRIFT WAVE 



A. The action functional 



In order to adimensionalize the equation (|l|) wc introduce the space and time scales t — > SI and y — > psy and the 
equation becomes 

^'-^-^Tt^\Ws)^y~\Ws)''^y^' ^''^ 

For simpHcity of notation we keep the symbol Vd for the adimensional velocity (^sfp") ■ "^^^ equation does not change 
of form but now all variables are adequately normalized and the action 

S = J dydt (16) 

is also adimensional. 

We have to calculate explicitely the scalar function x- Based on the extended knowledge developed in field theory it 
seems reasonable to assume that this function represents the generalization of the functions which have the opposite 
evolution compared to f. if ip evolves toward infinite time, then x comes from infinite time toward the initial time. 
If if diffuses then x anti-diffuses (see Ref. |3^). The general characteristics of this behaviour suggest to represent 
X as the object with the opposite topology than ip. If ip has a certain topological class, then x has the opposite 
topological class. If p is an instanton then x is an anti-instanton. In our case: if ip is the vortex solution, then x must 
the "anti-vortex" solution, with everywhere opposite vorticity compared to ip. In our case of a single vortex, x must 
simply be a negative vortex. 

In general terms, the direct {i.e. the vortex + random drift waves) solution p arises from an initial perturbation 
which evolving in time breaks into several distinct vortices (solitons) and a tail of drift waves, as shown by the Inverse 
Scattering Method. The functionally conjugated ("regressive") function x is at < = oo a collection of vortices and 
drift wave turbulence which evolving backward in time, toward t = 0, coalesce and build up into a single perturbation, 
the same as the initial condition of ip. We can restrict our analysis to the time range where the two functions has 
similar patterns (but opposite) which simply means to chose the time interval far from the initial and asymptotic 
limits. As shown by analytical and numerical studies, the vortices (positive and negative) are robust patterns and the 
time evolution simply consists of translations without decay. In conclusion we can take for the time range far from 
the boundaries t — and t = T: 

X = -P (17) 

To see this more clearly, we write down the action and then the Euler-Lagrange equations, with the current J 
included. 

Sj[x,v] = [ dy [ dtLj,^ (18) 
Jo Jo 



with the notation 



dip dip dip 
dt oy dy 



+ J^ip + J^x (19) 



When performing integrations by parts the boundary conditions of the two functions prevents us from taking the 
integrals of exact differentials as vanishing, but this just produces terms which do not contribute to the determination 
of the solution of extremum. We shall first change the Eq.(|l9|) such as to obtain by functional extremization an 
(Euler-Lagrange) equation for the function x'- 

r(i) i,dp dip ^'^,7,7 /on\ 

Lj,!p = - [V^xJ ^ + VdXg^ - ^dXV-g^ + + JxX (20) 

Now we write the condition of extremum for the action functional and obtain the Euler-Lagrange equation 

d 5L% ^ d SL% ^ d 5L% 5L% _ ^ ^^^^ 



dt 
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This equation can be written 



An equivalent form of the action is 



Sj[x,^] ^ I dy I dtLj,^ (23) 
Jo Jo 



with 



Lj^x = - fl^ + (V±</?j -K7 - Vd^^. + ^'d^ir: + Jv'P + -^xX (24) 



The equation Euler-Lagrange for the function x is obtained from the extremum condition on the functional Eq. (|2 



'^"^(n) 5x 

This equation reproduces the nonlinear vortex equation with an inhomogeneous term 



= (25) 



Comparing the homogeneous equations ( p2| ) (with J;^ 0) and ( p6[ ) (with = 0) we see that 

X = (27) 



is indeed the solution of the homogeneous equation ( |22| ) i.e. the negative vortex is the solution for x- 

We must remember that the "external" currents are arbitrary and later, after functional differentation, they will 

be taken zero. This allows us to start from the configurations given by the homogeneous equations and Eq.(p7|) and 

study the small changes using perturbative methods developed in the framework of the Inverse Scattering Transform. 

We will only use the current which will be denoted J and already take J^ — Q. 
The final form of the action which will be used later in this work is 

^ i-T r 2\ , dip dip 



Sj [x, ^] = / dy I dt\x{i- vi) + vdx^ - "'^^'^ + r (28) 



B. The condition of extremum of the action functional 



The Euler - Lagrange equations for the two functions x ^'^d ip are obtained from the first functional derivative of 
the action Sj: dSj/dx = and dSj/6ip — 0. The first equation (which is the original equation) has the solution 
(^. It does not depend on the current J (since the corresponding current has been taken zero). However, for 
uniformity of notation we shall write ipjs 

'PJs{x,y,t) = ips{x,y,t) (29) 
The second Euler-Lagrange equation is the equation for x, with the inhomogeneous term given by the current J: 

(,_Vi)|..|-..|.. ,30, 

The solution is: 

Xjs {x, y,t) = - ifs [x, y, t) + xj {x, y, t) (31) 

where — ips {x,y,t) represents the "free" solution of the variational equation, i.e. the negative vortex (anti-soliton) 
and Xj {x,y,t) is the small modification induced by an inhomogeneous small term, J{x,y,t). Since the function 
Xj (x, y, t) is the perturbation of the negative-vortex solution we will use the equation ( p^ ) but with the opposite 
current {i.e. —J instead of J), as (pO|) requires. 
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C. Second order functional expansion and the eigenvalue problem for the calculation of the Determinant 



Now we shall expand the action Sj [ip] to second order around the saddle-point solution. Write 



ip = (pjs+ Sip 
X = XJs+ 6x 



(32) 



where the function (5ip, Sx) is a small difference from the extremum solution. The expanded form of the action will be 
written: 



Sj [x, = S.J [<pjs,X.Js] + 7:1 T^. 



\ S(pSx 



SipSx 



V.Jb.XJb 



where obviously the absence of the linear term is due to the fact that {ipjs, XJs) is the solution at the extremum and 



Sj [vjs,XJs 



dy / dt 



XJs 



dip. 



Js 



dt 



f 2 \ dfjs 



(33) 



dipjs dipjs , 

+VdXJs—^ VdXJsP>Js—^ h Jipjs 



dy 



dy 



Few manipulations are necessary to make the second functional variation of Sj symmetric in 5ip and 5x- Again 
this will imply boundary terms, but these are now zero since the variations Sip and 5x vanishes at the limits of the 
space-time domain, by definition. The transformations are simply integrations by parts and give 



1 



■5X 



S^Sj 



\ SipSx 



VJs,XJs 



Sip^^ Sip Sx 



7 —a — (3 
a-P 



Sip 
Sx 



where 



S = 1 - V 



d 



d 



f3 = -vd 



dt 

dipjs 
dy 



dy '"'^ V dy 



( dip 



Js 



(34) 



(35) 



^ o d 
7 = -^VdXJs^ 
dy 

In the generating functional of the correlations, the expansion gives, after performing the Gaussian integral: 

-1/2 



1 



Zj = exp(i5j)-(2^) 



\n/2 



det 



(S^Sj 



\ SipSx 



As stated before, the det will be calculated as the product of the eigenvalues A„ 



det 



( S'Sj 
\ SipSx 



We must find the eigenvalues of the differential operator appearing in Eq.(p4[) 
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a- 13 



which gives the following equation 



(36) 



(37) 



(38) 



(39) 
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The functions Sip (y. f.) and Sx (?/■ t) represent the differences between tlie solutions at extremum (solitons) and other 
functions which are in a neighbourhood (in the function space) of the sohtons. According to the discussion above, 
the functions which are "close" to the solitons, for the Flier 1- Pet viashvilli equation are drift waves. For this reason 
the operator which represents the dispersion {i.e. V^) will be replaced with its simplest form , — fc^ for these waves, 
with k± representing an average normalized wavenumber for the pure drift turbulence. However, the operator will be 
retained when applied on the functions related to solitons, since these solutions owe their existence to the balance of 
nonlinearity and dispersion. The following detailed expressions are obtained for the operators involved in this equation: 



a 



d 



d 



dt 



dy 



-Vd 



(1-vi) 



dt 



d_ 

dy 



+Vd (1 - 'PJs) (l + fcj. 
-vl (1 - <fJs) 



92 



dip 



Js 



dy 



(40) 



dtdy 
d(pjs\ d 



dy J dy 



Vd 



(1 - fJs) (1 



dydt 



- 1 / d(pj. 



dt 



(1-vi) 



dtdy 



Js 



dy"^ 



dy 



(1 - V>Js) 



dy 



(41) 
(42) 



d<fjs\ d 



dy J dy 



d(pj, 
dy 



(43) 



The square brakcts are used to underline that the differential operators are not acting outside and the only operation 
is multiplication. We use the equation verified by ipjg to make the following replacement 



The equation becomes 



(1-vi) 



dt 



-Vd (1 - ifjs) 



dip 



Js 



dy 



- I (1 + fcl)' ^ + 2 (1 + fcl) Vd (1 - fjs) 



dydt 



(44) 



(45) 



,2/1 ^2 92 82/ d<pjs 



d 



Tn 



+Xn[-2VdXJsg^]^^ 



We now take into account the propagating nature of the drift waves and make the change of variables t ^ t 
and y ^ y — Vdt i.e. we change to the system of reference moving with the diamagnetic velocity. We simplify the 
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equation assuming that the most important space-time variation is wave-Uke and replace ^ = —Vd-§jj- By this change 
of variables the soliton will not be at rest in the new reference system, but it will move very slowly since we have 
assumed that Vd- We make another approximation by neglecting the slow motion of the soliton. This restricts us 
to the wavenumber spectrum but considerably simplifies the calculations. The space variable which will be denoted 
again y measures the space from the fixed center of the soliton, in the moving system. The difference between the KdV 
soliton, which is one-dimensional and depends exclusively on y and the vortex which is a two-dimensional structure 
will be considered in the simplest form as described by the estimation of Meiss and Horton for the x - extension of 
the vortex. For convenience we suppress the index n and replace ipi^ by q. 



1 + k^] Vd- I'd (1 - tpjs) 
d 



2 



(46) 



{2XvdXJs) 



dy 
( dip 



Js 







We have a suggestive confirmation that the generating function Zj (via the action Sj) potentially contains con- 
figurations of the system consisting of simple drift waves. A perturbation consisting of drift waves and propagating 
with the diamagnetic velocity Vd is an approximate solution of the original equation for small amplitude (i.e. small 
nonlinearity) . Due to its particular structure, the Martin-Siggia-Rosc action functional is exactly zero when calculated 
with the exact solution, in the absence of any external current J. The action expanded to the second order then gives, 
for no vortex {cpjs = 0, xjs = 0) 



A 



k^,Vd, 



q = 



which implies periodic oscillations in the space variable y with (recall that everything is adimensional) 

X = kyVd (k]_ 



Returning to the equation (^ ) we write it in the following form 



where 



A = 



B 



2A X 



Js 



(t? 



^ 4 



(47) 



(48) 



(49) 



(50) 



Now we make the standard transformation of the unknown function 

g = wexp (-]^ [ A{y')dy' 



and obtain 



w"+(B-4--^|w; = 



(51) 



(52) 



where prime means derivation with respect to y. After replacing the two extremum solutions ipjs and xJs from 
equations (29) and ( ^ this equation is written in the following form, to exhibit the dependence on A: 
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w" + {X^ti + Xt2+t3)w = (53) 

with the notations 
and 



4 /i2 V 9y 



c = fcl (57) 

h — c + ^Ps 

The functions (y) are represented for i = 1, 2 in Figure 2 and 3. The function 

U {\;y) = \^ti+\t2+h (58) 

has singularities at the points where h vanishes. We introduce the notation yh for the location of the singularities, 
taking into account the symmetry around y = 0, the centre of the soliton 

h {±Vh) - 0. (59) 

Since the soliton is very localized, the function U has very fast variations close to the singularities. The slow variation 
of the function U {X;y) over most of the space interval {—L/2,+L/2) becomes very fast due to the growth of the 
absolute values of ti, t2 and near ±y/i, on spatial intervals having an extension of the order of the spatial unit, 
i.e. ps in physical terms. Since the physical model leading to our original equation cannot accurately describe the 
physical processes at such scales, we shall adopt the simplest approximation of U , assuming that it reaches infinite 
absolute value at points which are located whithin a distance of ps of the actual positions of the singularities, -iiyh- We 
have checked that the exact position of the assumed infinite value of U has no significant impact on the final results, 
which can be explained by observing that ti,2,3 will be integrated on. The total space interval is now divided into 
three domains: {—L/2,—yh) (external left), {—yh,yh) (internal) and {yh,L/2) (external right). Here "internal" and 
"external" refer to the region approximatly occupied by the soliton. The form of the function U imposes the function 
w to vanish at the limits of these domains. In a more general perspective, the fact that w behaves independently on 
each domain has a consequence with statistical mechanics interpretation: the generating functional (similar to any 
partition function) is obtained by integrating over the full space of the system's physical configurations and behaves 
multiplicatively for any splitting of the whole function space into disjoint subspaces. In particular the functional 
integration over the space of functions Sip and 6x actually consists of three functional integrations over the disjoint 
function subspaces corresponding to the three spatial domains. The fact that our physical model is restricted to 
spatial scales larger than ps necessarly has an impact on the maximum number of eigenvalues A„ that should be 
retained in the infinite product giving the determinant, but we shall not need to use this limitation. 

For absolute values of the parameter A greater than unity (which will be confirmed a posteriori, by the expressions 
( |6^ ) and (^l|) below), the three terms in the expression of U have very different contributions. The terms ^3 is 
practically negligible, and the term with ti is always much greater than t2 in absolute value. In the following we 
consider separately the three domains. 

On the "external left" domain, the function ti is positive. If we fix at zero the amplitude and the phase of w at the 
limit —L/2 the condition that the solution vanishes at —yh gives, for A real, 

^'^ dy' (A^ ti+Xt2 + t^y^^ = 27rn (60) 

L/2 

In the integrand, the first term is factorized and, taking into account the relative magnitude of the terms, we expand 
the square root and obtain 
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i.e. 



Xlai + A + ^ = 27m (61) 



= ^ f 1 - MM , (62) 



ai 



where 



ai 



dy' VhW) (63) 

L/2 



L/2 Vil (y') 



Ti = / /-^^ (65) 



V2 VW) 

and 7i has been neglected. We note that /3i is positive. 

On the "external right" domain the function ti is positive but t2 is negative. The condition on the phase is 



/ dy' {\^ti+\t2+uY''^ = 2-Kn' (66) 

J Vh 



and introduce similar notations 



a2 



L/2 



dy'Vti {y') = ai (67) 



Vh 



The equation then becomes 



/32= r'dy'-^Hi^ = (68) 
h jy' 



A;;,a2+/32 + ^ (70) 



r'\y'-h^ (69) 



or 



Q!2 V / 

The infinite product of eigenvalues gives, for the "external" region |S: 



n n' n ^ ' n \ / 

_ sin(/3i/2) -|-r /27rn 
" /?i/2 y 

In the "internal" region, the function ti is negative. The relations between the magnitudes of the absolute values 
of the functions ii, t2 and are preserved. Then A will be complex. Due to the anti-symmetry of the function t2 we 
can suppose that the unknown function w takes zero value at y = 0. We introduce the notations 
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Jo 



(73) 



= r\y' 



which are rea/ numbers. The condition 



(2/0 



A>e + /3c + = 27rm 



gives (after neglecting 7c) for the complex number X\: 



PI 



1/2 



At, = (27rn) 1 

\ (27r) 



The infinite product of these eigenvalues is 

n^n = n"c"'(27rn)exp 



-i arctan 



exp 



2?™ 



—i arctan 



n ' 



/27rn 



/3cV(27r)^ 



(74) 



(75) 



(76) 



(77) 



1/2 



(78) 



The number /3c is smaller than unity and for large n the argument of the exponential will be more and more close to 
— i7r/2. We make the approximation that the exponential can be replaced with —i. Then we obtain 



sinh(/3c/2)]^/^j-|- (-i)27rn 



/3c/2 



an 



(79) 



On the "external" regions the functions ti, t2 arc not symmetrical around the centre y = since the perturbed 
soliton develops a "tail" which is not symmetrical. However we take this perturbation to be small and assume the 
same absolute value for the function f3i on both external domains. 

Wo remark that wc remain with two quantities in which all the functional dopcncence on the current J is packed: 
for "exterior" /?i (hereafter denoted a) and for "interior" (3c (hereafter denoted /3). 



const ex.p{iSj) 



det 



-1/2 



(80) 



r /3/2 1 


1/4 


a/2 


1/2 


_sinh(/?/2)_ 




sin (cr/2) 





where 



const = 



{-i)ac 
2'Kn 



1/2 



n 



(81) 



will disappear after the normalizations required by the calculation of the correlations (see below). 



V. CALCULATION OF THE CORRELATIONS 

The two-point correlation can be obtained by a double functional differentiation at the external current J. 

5^Zj 



(^(t/2)<^(yi)) = ^7^ 



iSJ{y2)iSJ {yi) 



j=o 
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The main achivement of this approach is that it provides the exphcit expression of the generating functional. We 
introduce the notations 



A = A[J] 



/3/2 
sinh (/3/2) 



1/4 



(82) 



B = B[J] = 



ct/2 



sin (ct/2) 



1/2 



(83) 



and drop the factor const] actuaUy the latter depends on ai and Uc and thus on the current J and contributes to the 
functional derivatives. However, taking a formal limit N to the number of factors in (|^) we find that the functional 
derivatives of ai and ac give additive terms which vanish in the limit TV oo. Then we drop const since it disappears 
after dividing to Zj and taking J = 0. In this way (|80| ) becomes 



Zj = exp(iS'j) AB 



(84) 



We calculate the functional derivatives. 



SZj 



iSJ (yi) 



SS., 



1 SA 



1 6B 



5J{yi) AiSJiyi) B iSJ (yi) 



exp(i5j) AB 



We will also need the functional derivative at J (1/2)7 with a similar expression. The second derivative: 



S^Zj 



i6J{y2) 16 J (yi) 



SSj SS, 



SJ{y2)SJiyi) i5J{y2)5J{y^) 



(85) 



(86) 



1 5A SSj 
'AiSJ{y2)SJ{yi) " 

1 6A SSj 
'AiSJ{yi)SJ{y2) ' 

1 SA I SB 
'AiSJ{yi)BiSJ{y2) 

1 S^A 



1 SB SSj 
BiSJ{y2)SJ (yi) 
1 SB SSj 
BiSJ{yi)SJ iy2) 

I SA 1 SB 
'^AiSJ {y2)BiSJ{yi) 
S^B 



Ai5J{y2)iSJ{yi) B iSJ {y2)iSJ {yi) 

The detailed expressions of these terms are given in the Appendix. The terms are calculated numerically using the 
detailed expressions of Lps, ^§^, XJ 

The first term reproduces the self-correlation of the soliton and represents the connection with the results of Ref. 
I p^ , with our particular simplifications: single soliton and fixed (non-random) position of its centre. As can easily 
be seen, the first order functional derivatives of Sj to the current J reduce to the function ips calculated in the 
corresponding points. The term with the double functional derivative of the action represents the contribution to the 
self-correlation of the soliton due to a statistical ensemble of initial conditions, without drift waves. All mixed terms 
{i.e. containing both the action and one of the factors A or B) represent interaction between the perturbed soliton 
and the drift waves. The terms containing exclusively the factors A and/or B refers to the drift waves in the presence 
of the perturbed soliton. 



VI. DISCUSSION AND CONCLUSIONS 



The formulas obtained by functional differentiation of the generating functional are complicated and a numerical 
calculation is necessary. We chose a particular value of the soliton velocity (which also fixes its amplitude): u = 1.725 Vd 
and let the variables yi and ?/2 sample the one-dimensional volume of length L = 0.2 m. The physical parameters are 
chosen such that ps ~ 10~^m and Vd ~ 571 m/s. We recall that there are two particular symmetry limitations of our 
calculation. (1) The soliton centre is assumed fixed (at y = 0) , especially for avoiding the complicated problem of 
the zero modes. (2) Due to the asymmetry of the perturbed soliton tail the terms which results from the functional 
differentiation are also asymmetric. These are only limitations of our calculation and in no way reflect the reality of a 
isotropic motion of many solitons in a real turbulent plasma. In order to see to what extent our result can be useful 
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for understatnding the (much more complicated) real situation we will symmetrize these terms in the unique mode 
which is accessible to our one-dimensional calculation, i.e. take into account the mixing of perturbed solitons moving 
in the two directions on the line. 

The amplitude of the modifications of the soliton depends on a parameter which is the average time of interaction 
with the perturbation. This average time is comparable with the time required to cross L at a speed of Vd and is 
limited since the growth of the perturbation cannot exceed the soliton itself. 

The figures are conventional representations of functions of two variables (j/1,2/2); they do not correspond to a 
two-dimensional geometry. For this reason it is not expected to have circular symmetry. The contributions to the 
correlation from the last two factors in Eq.(|8^) have amplitudes similar or less by a factor of few units, compared to 
the pure soliton. The factors coming from "internal" part are peaked and localized on the soliton extension while the 
"external" part gives terms oscillating on (2/1,2/2)- In wavenumber space, there are contributions to both low-fc and 
high-fc regions. The spectrum of an unperturbed soliton is smooth and monotonously decreasing from the peak value 
at k = 0. Fig. 5 shows much more structure. In the low-fc part there are many local peaks, an effective manifestation 
of the periodic character of the terms (as shown by (72) ). This arises from the discrete nature of the eigenvalues, 
which is induced by the second order differential operator and the vanishing of the eigenmodes at the positions of the 
singularities « ±yh. The singularities are generated by the vanishing of the norm of the operator a, which makes 
ambiguous the assumption of propagating wave character, dt — —Vddy. The large-fc part mainly reflects the structure 
of the small-scale shape perturbation of the soliton, comming from /3-related terms. Fig. 6 is an (fc, w) spectrum 
obtained from lo — ku = Q and repeating the calculations for various soliton velocities Umax > u > Vd- Although we 
cannot afford high Umax since the expressions of ii,2,3 {y) depend on the assumption u Ud , we remark local peaks 
in contrast to the "pure soliton" result of Ref. H. 

For simplicity we have assumed a single soliton. However the calculation can be readily extended to the multi- 
soliton case, considering instead of ( p9| ) and ( ^l| ) sums over many individual soliton solutions with different velocities 
and positions of the centres. These sums replace the functions ipjs and XJs in the expressions of the operators 5?, (3 
and 7. If the velocities are all greater but not too different of Vd the change of variables to the referential moving 
with Vd (described in the paragraph below Eq.(^) ) will leave a very slow time variation which eventually may be 
treated perturbatively. Many solitons will also generate many singularities arising from the vanishing of the function 
/i, and this will factorize the space of functions and correspondingly the generating functional. It will become however 
possible to consider random positions and random velocities and average them with distribution functions for the 
Gibbs ensemble, like in |p^ . This is very simple with the first term of (|86|), which should be compared directly with 
Ref. [ p^ , but technically very difficult with the terms involving functional derivatives of A and/or B. 

The first results suggests that the non-gaussianity at the plasma edge can be explained by the presence of coherent 
structures. The contribution of avalanches to the deviation from the gaussian statistics cannot be excluded but, as 
shown for self-organized systems | ]36| , they have a scaling which should be easy recognized, at least in frequency 
domain. 

In conclusion we have developed an approach which allows us to calculate the statistical properties of a coherent 
structure in a turbulent background. Compared to the standard renormalization, this approach is at the opposite limit 
in what concerns the relation "coherent structure / wave turbulence" , highlightning the coherent structure. However 
it offers comparatively greater possibilities for the extension of this studies to the more realistic problem of cascading 
wave turbulence mixed with rising and decaying coherent structures. 
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APPENDIX A: EXPLICIT EXPRESSIONS FOR THE FUNCTIONAL DERIVATIVES 

We shall first concentrate on the derivatives of the two factors A and B. 



5B 



6J{yi) 5J{yi) 

_ 1 
~ 4 



0-/2 



■ 1/2 
sin (cr/2)_ 

fj cos {a/2) 



sin (a/2) 2sin2(cr/2) 



a/2 



sin {a/2) 



(Al) 



-1-1/2 



5a 



SJ iyi) 



and 
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SJ{y2)SJ " SJ{y2)SJ{yi) 



a/2 
sin (a/2) 



1/2 



1 

~8 
1 

16 



a/2 


1/2 ^ 


_sin(cr/2)_ 




0-/2 


-3/2 r 


sin (a/2)_ 





+- 



For the exterior domains, 



a/2 


-1/2 


1 


sin ((t/2) 




sin {a/2) 



1 + cos^ ((t/2) 
sin^ {a/2) 

1 cr cos {a/2) 

sin(<T/2) ^ 2sin2 {a/2) 

a cos {a/2) 



} (<5J(2/2)) (5J(2/i) 
\SJ{y2)SJ{yi)J 



with 



a = ao + aji + aj2 



CTO 



1 /■"'"' 

= 2/. *' 



-Z,/2 



f-1 



dy ) (/i2_ ^2)1/2^ 



(A2) 



(A3) 



(A4) 



L/2 



dy' 



dy J /i(/i2- ^2)1/2 V K'-f . 



(A5) 



crj2 



1 Z""^'' 



-t/2 



1 /5x7 

(/l2- ^2)1/2 1^ dy 



We have the following connected expressions: 

5a 



5a. 



71 



5a 



72 



(5J(yi) W(2/i) 5J{yi) 



(A6) 



(A7) 



and: 



<55 



ji 



<525ji 



^J(t/2)<5J(2/i) 



SJ{yi) 2j_L/2 ' V dy J h{h?-^lf'^ 

(-1) 



2- 



5?, 



J2 



(2c - h) 



SJ {yi) 2 J_L/2 (/i2 _ <p2)i/2 6J {yi) V dy 



5^a 



5^ 



CTji 



''■,72 



5J{y2)5J{yi) dJ{y2)5J{yi) dJ{y2)5J{yi) 

7) = 2L2''[-W) 



mm) 



_ {2c-h) \ ( 5^x7* \ 
h{h^~^lf'^\ J\5J{y2)5J{yi)J 



(A8) 



(A9) 



(AlO) 



(All) 



dy' 



(-1) 



6J{y2)6J{yi) 2 J_j^/2 {h^ - vlf'^ 5J{y2)5J{yi) \ dy 



(A12) 



For the "interior" region, the derivatives of A, (which are strightforward) will require the calculation of the deriva- 
tives of p. 
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1 rVh 
/3= o 



J_ ( 2c-h I 2_ 1 f d!£^\ ~int ]_ 1 rfx!" 



( 1 ^l~h^ \^^'^ ( 



_2Vs_~int\ 



1/2 



The function x'"* and its derivative are present in the expression of /3: 

P = P0+PjI+Pj2 



1 f^'^ 

1^0 = t; dy 



difs \ 2c — h 



1 f^'^ 



dip, 

dy J h {ifl - ft2) 



1/2 



ips (2c - h) , 



1 fy^^ 

(ij2 = -,l dy 



dxT 

(¥^2 _ /i2)l/2 y 



and the derivatives at J arc easily calculated, as for a. 

The formulas above need to specify the expression of the functions xT* > ^^y ^^'^ '^^ their functional derivatives. 
We use the results of the analysis carried out by Karpman. 
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Figure Captions 



Fig.l The form of the soliton fs{y) for the velocity u = 1.725 u^- 

Fig.2 The function ti {y) for u = 1.725 w^- 

Fig. 3 The function t2 (y) of the Eq.(|5^) for the same u. 

Fig. 4 The perturbation to the correlation in physical space. 

Fig. 5 Contour plot of the spectrum of the vortex perturbed by the turbulent drift 
Fig. 6 The contour plot of the frequency-wavenumber spectrum, with ui — ku = 0. 
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FIG. 1. Fig.l Variation of the form of the soUton ifs{y) with the velocity, u. 
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FIG. 2. Fig. 2 The function ti (y) for a particular soliton velocity, u = 1.725 1;^. 
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FIG. 3. Fig. 3 The function t2 (y) of the Eq.(B3) for the same u. 



20 




FIG. 4. Fig.4 The perturbation to the correlation in physical space. 
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FIG. 5. Fig.5 Contour plot of the spectrum of the vortex perturbed by the turbulent drift waves. 
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FIG. 6. Fig.6 The contour plot of the frequency-wavenumber spectrum, with uj — ku = 0. 
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